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In this article we discuss the realization of atomic GHZ states involving three-level atoms in a 
cascade and in a lambda configuration and we show explicitly how to use this state to perform 
the GHZ test in which it is possible to decide between local realism theories and quantum 
mechanics. The experimental realizations proposed makes use the interaction of Rydberg atoms 
with a cavity prepared in a state which is a superposition of zero and one Fock states. 
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X : 1 INTRODUCTION 

a' 



The superposition principle, entanglement and its consequence non-locality makes quantum mechan- 
ics fascinating and, at the time, intriguing and not intuitive according to our common sense which 
is based essentially on classical mechanics. Quantum mechanics has given rise to many questions of 
philosophical nature. For instance, the concept of reality PjJ El HI E] has new interpretations under 
the light of quantum physics. The Consequences of entanglement and non-locality were noticed by 
Einstein, Podolsky and Rosen who proposed a gedanken experiment, the EPR experiment in a very 
interesting article [HI 121 El HI Ej to show that quantum theory was not a complete theory. In their 
article they describe an experimental arrangement involving correlated pair of particles where these 
particles interact and then are separated. .Therefore, measurements made on one particle can be 
used via correlation to generate predictions about the other particle. EPR state that "If, without 
in any way disturbing a system we can predict with certainty the value of a physical quantity, then 
there exists an element of physical reality corresponding to this physical quantity" and as a necessary 
condition for the completeness of the theory, EPR state that "every element of the physical reality 
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must have a counter part in the physical theory" . EPR noticed that the fact that quantum mechan- 
ics admits no dispersion- free states, does not tell us wether the theory is complete or not. In their 
words: "From [the dispersion principle] it follows (1) the quantum mechanics description of reality 
given by the wave function is not complete or (2) when operators corresponding to two physical 
quantities do not commute the two quantities cannot have simultaneous reality" . Shortly after the 
EPR paper Bohr published a reply in which he defends the completeness of the quantum-mechanical 
description of nature which could not be refuted properly by EPR and which came to be known as 
the Copenhagen interpretation jS]. 

Hidden variable theories [21 E] were proposed in which one would recover local realism and one 
of its proponents was Bohm [OlEllSj- Bell developed a clever theoretical tool which could be used 
to perform experiments based on it which could decide between local realism theories and quantum 
mechanics, the Bell's inequality [TUl 13 HH Hj- In the 1980s it was realized experiments by Aspect 
and collaborators fTJ [7j based on the Bell's inequality which although strongly favored quantum 
mechanics, there remained some possibility for which a local reality view could still be maintained. 
We should point out that the main difficulty in the test of the Bell's inequality is related to the fact 
that there is no perfect particle detectors. For a review of Bell's inequalities and some of its variants 
see P2J. 

In 1989, Greenberger, Home and Zeilinger found out an ingenious theoretical tool to test quantum 
mechanics confronted with local hidden variable theories, the Bell theorem without inequalities [Till 
HUH]], which can demonstrate the spookiness of quantum mechanics even more dramatically than the 
Bell's analysis. In the GHZ test the decision between quantum mechanics and local realism theories 
is a binary simple one showing the power of the GHZ formalism. There are several proposals of 
preparation of GHZ states. For instance, proposals involving cavity QED are presented in |T71 EH 
IT8*1 For a concrete experimental realization of GHZ states see |19| ITo] . In reference we 
have proposed a preparation of atomic GHZ states and a GHZ test similar the one discussed in this 
article. However, in reference ^7] we have used three-level lambda atoms in which the lower states 
are degenerate states interacting with a cavity prepared in a coherent state. In the present article 
we assume that the atoms are Rydberg atoms of relatively long radiative lifetimes [20J and that the 
superconducting microwave cavities [2U [22] are perfect cavities, that is, we neglect effects due to 
decoherence. We are going to consider cavities which are prepared in a superposition of zero and one 
Fock states. 



2 ATOMS IN A CASCADE CONFIGURATION 



We start showing how to prepare a cavity C in the state 

(|0> + |1» 



V2 



(2.1) 



In order to prepare this state, we send a two-level atom AO, with |/ ) and |eo) being the lower and 
upper level respectively, through a Ramsey cavity RO (see (|A.72jl ) in the lower state | fo) where the 
atomic states are rotated according to 



Ro = 7^ 



1 i 
i 1 



that is, 



fo) -> -j=(i I e )+ | fo)), 



(2.2) 
(2.3) 
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and through C, resonant with the cavity. Under the Jaynes-Cummings dynamics [23] (see (jA.67|) . 
for A = 0). For gt = 7r/2 we see that the state |/o)|0) does not evolve, however, the state |eo)|0) 
evolves to — i|/o)|l). Then, for the cavity initially in the vacuum state |0), we have 



(|/o> + i|e >) 
V2 



10) — |/o> 



(|0> + |1» 
V2 



\foM 



c- 



(2.4) 



Now let us consider a three-level cascade atom Ak with | e^), | fk) and | g^) being the upper, 
intermediate and lower atomic state respectively (see Fig. 1). We assume that the transition | /&) ^| 
Ck) is far enough from resonance with the cavity central frequency such that only virtual transitions 
occur between these states (only these states interact with field in cavity C). In addition we assume 
that the transition | e^) ^| gk) is highly detuned from the cavity frequency so that there will be no 
coupling with the cavity field. Here we are going to consider the effect of the atom-field interaction 
taking into account only levels | fk) and | gk)- We do not consider level | e^) since it will not play 
any role in our scheme. Therefore, we have effectively a two-level system involving states | fk) and 
\gk). Considering levels | fk) and | gk), we can write an effective time evolution operator 



u k {t) 



fk)(fk I +\9k){gk 



(2.5) 



(see (IA.70j) ) where the second term above was put by hand just in order to take into account the 
effect of level | gk) and where <p = g 2 r / A, g is the coupling constant, A = u e — ujf — u is the 
detuning where u e and Uf are the frequencies of the upper and intermediate levels respectively and 
to is the cavity field frequency and r is the atom-field interaction time. Let us take (p = ir. Now, let 
us assume that we let atom Al to interact with cavity C prepared in the state (j2.1J) . Let us assume 
that atom Al is sent to a Ramsey cavity Rl in the lower state | gi) where it is prepared cprepared 
in a coherent superposition according to the rotation matrix 



and we have 



1 

1 



1 1 
-1 1 



1>)Ai=j={\fi)+\gi)). 



(2.6) 



(2.7) 



Taking into account ()2.5j) . after atom Al has passed through the cavity prepared in state (j2.1)l . we 
get 

I ^)ai-c = A>+ I 9i))\0) + (- | A>+ I <7i»|l>], (2.8) 

Now, if atom Al enters a second Ramsey cavity R2 where the atomic states are rotated according 
to the rotation matrix 



we have 



> 



and, therefore, 



1 


i i ' 




= V2 


-i i 




1 A>- 


t- 1 9i)) - 


I/O, 


1 h)- 


f 1 9i)) - 


1 9i), 


= T2 [ 


/i>|0>+| 


0l>|1>] 



(2.9) 



(2.10) 



(2.11) 
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Now, let us prepare a two-level atom A2 in the Ramsey cavity R3. If atom A2 is initially in the state 
| #2), according to the rotation matrix 



we have 



j_ r 1 1 

y/2 [ "I 1 
V^)a 2 = -^=(| / 2 )+ I <fe)), 



(2.12) 



(2.13) 



and let us send this atom through cavity C, assuming that for atom A2, as above for atom Al, 
the transition | f'2) ^| ^2) is highly detuned from the cavity central frequency. Taking into account 
(|2.5|) . after the atom has passed through the cavity we get 



A1-A2-C 



;\1 fi)(\ f*)+ \ 9*W)+ \ 9i)(- \ f*)+ \ 



(2.14) 



Then, atom A2 enters a Ramsey cavity R4 where the atomic states are rotated according to the 
rotation matrix 



Ra — — F= 



1 

V2 



1 1 
-1 1 



(2.15) 



that is, 



and we get 



A1-A2-C : 



1 / 2 >+ 


1 92)) - 


♦1 h), 




1 /2> + 


1 <72» - 


*l 92}, 


(2.16) 



1 



T2 [ 



fi) |/ 2 >|0>+|pi> I fla>|l>], 



(2.17) 



which is a GHZ state. If we had started with the cavity in the state (|0) — \ l))/y/2 we would get the 
GHZ state 

\4>,-)M-A2-c=A=[\fi) \f2)\0)-\ 9l ) \g 2 )\l)}. (2.18) 

Let us now first discuss a summary of the GHZ test prescription. We will follow closely the 
discussion presented in the very clear article by Mermin jT3]. First we define the atomic operators 



A 
B 
C 
D 



12 3 
(J x (T y (J y^ 

^-1^-2 3 
12 3 
^1^2 3 



(2.19) 



where 



/fcXfl 1 * I + I 9k) (fk \, 
~K\ fk)(9k I - I 9k) (fk |), 
0)(1| + |1)(0|, 
0)(1 I - I 1)(0 |), 



(2.20) 
(2.21) 
(2.22) 



(k = 1, 2 ). It is easy to show that the commutators 

[A B) = [A, C] = [B, C) = 0, (2.23) 

and that 

A | V, ±)a1-A2-C = B | V, ±}a1-A2-C = C | V, ±)a1-A2-C = =Fl I i> ■> ±)ai-A2~c, (2.24) 

and 

D | V, ±)ai-A2-c = ±1 I 4, ±)ai-A2-c- (2.25) 

If we assume that there are six elements of reality m k x and m k y (k = 1,2 and 3) each having value 
+1 or —1 waiting to be revealed according to a local realism theory, then we can write 





= m\m\m\ 


= Tl, 


b± 


= m\m\m\ 


= Tl, 


c± 


= m\m\m\ 


= Tl, 



(2.26) 
and we have 

d± = a±b±c± = m l x m 2 x ml = =f1, (2.27) 

where the upper and lower subindexes refer to the GHZ state (j2.17j) and (j2.18j) respectively and we 
have used {rn k ) 2 = 1. So, the existence of elements of reality implies that if we measure the value of 
the observables a x (k = 1,2 and 3) (that is the elements of reality associated with them) in the state 
| ±)ai-A2-c the product of the three resulting values must be d± = =Fl. But according to ()2.25|) 
the eigenvalue of the operator D applied the state | if), ±)ai-A2-c is ±1. Therefore, measuring this 
eigenvalue we can decide between theories based on local realism and quantum mechanics. 

Now, let us see how we proceed to perform the GHZ test (we follow closely the scheme presented 
in ^7j). We start letting atoms Ak to pass through the Ramsey zones Kk (k = 1,2 ) where the 
atomic states are rotated according to the rotation matrix 



K> = 75 



1 -1 
1 1 



(2.28) 



or 

K k = -^=(| f k )(f k | - | f k )(g k | + | g k ){f k | + | g k )(g k |). (2.29) 

These rotation matrixes are the key ingredient which allow us to perform the GHZ test. The method 
of the test we are going to describe is based on a gradual unraveling of the GHZ state being considered. 
The eigenvalues of the operators a x are 

|V&±) = ^(| A>±|0*»- (2.30) 

for k = 1,2. Let us take the state f)2. 17j) . Writing the states | f\) and | g\) in terms of the states 
()2.30|) for k — 1 and substituting in the GHZ state ()2.17|) we get 

I ^+)ax-A2-c = ~[|^,+>(l h) I 0)+ I 92) I 1)) + I A -XI h) I 0)- I g 2 ) I 1))]. (2.31) 
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Applying (|2.29j) in the state (|2.17|) for k — 1, we have 

K x | ^+)ai-A2-c = \[\fi)(\ / 2 ) I 0)- | g 2 ) | 1)) + 1^(1 f 2 ) | 0)+ | g 2 ) | 1})]. (2.32) 

Now, we compare ()2.32|) and (j2.31|) . We see that the rotation by K\ followed by the detection of |<?i) 
corresponds to the detection of the the state +) whose eigenvalue of o\ is +1. After we detect 
we get 

I ^)A2-c = ^=(1 / 2 > I 0)+ I g 2 ) | 1)). (2.33) 

As before, we rewrite the states | f 2 ) and | g 2 ) in terms of these states ()2.30j) for k = 2 and substitute 
in the above state and we have 

I ^)A2-c = \UI +>(| 0)+ | 1)) + |V£ ->(| 0)- | 1))]. (2.34) 

If we apply fl2~29l for k = 2 to the state (l2~33l we get 

K2 | ^)a2-c; = i[|/ 8 >(| 0)- | 1)) + |^ 2 )(| 0)+ | 1))]. (2.35) 

Again, if we compare ()2.35|) with ()2.34|) . we see that the rotation by K 2 followed by the detection 
of | g 2 ) corresponds to the detection of the the state \ipl, +) whose eigenvalue of is +1. After we 
detect \g 2 ), we get 

h/4+> = ^(|0>+|l», (2.36) 

which is the eigenvector of a z x with eigenvalue +1. We now send a two- level atom A3, resonant with 
the cavity, with \g 3 ) and I/3) being the lower and upper level respectively, through C. If A3 is sent 
in the lower state ^3), under the Jaynes-Cummings dynamics, for gt = tt/2, we know that the state 
|g 3 )|0) does not evolve, however, the state |g3)|l) evolves to — «|/ 3 )|0). Then we get 

lg , m+) ^M_m m (237) 

Now we let atom A3 to enter a Ramsey cavity K3 where the atomic states are rotated according to 
the rotation matrix 

" 1 i 



i 1 



(2.38) 



that is, 

I h)+ i<?3» H</3>, (2.39) 

and we see that the rotation by K 3 followed by the detection of \g 3 ) corresponds to the detection of 
the the state +) whose eigenvalue of a\ is +1. 
In the case of 

|1&-> = ^(|0>-|1» (2.40) 
which is the eigenvector of a\ with eigenvalue —1, we would get 

Ig3)fe-)^ te) ^ l/3>) |0) (2-41) 
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and after the rotation (J2.38|) we would have 



^(i|/3>+|fti»-i|/3>. (2-42) 

and we see that the rotation by K 3 followed by the detection of |/ 3 ) corresponds to the detection of 
the the state — ) whose eigenvalue of a^. is —1. 

We can repeat the above procedure and see that we have only four possibilities which are presented 
schematically below, where on the left, we present the possible sequences of atomic state rotations 
through K k and detections of | /_,•) or | gj) (J = 1,2) and | g 3 ) or | / 3 ) and on the right, we present 
the sequences of the corresponding states \if)%, ±) where k = 1, 2 and 3 which corresponds to the 
measurement of the eigenvalue of the operator D given by ()2.19j) . 





1 9i))(K 2 ,\ 


92))(K 3 ,\ 


9s)) «- 






(#i, 


1 9i))(K 2 ,\ 


h)){K 3 ,\ 


/s»<" 








1 fi))(K 2 ,\ 


h)){K*A 


9a)) <~ 




-M,-M,+>, 




1 fi))(K 2 ,\ 


92))(K 3 ,\ 


/s»<- 


-Mi, 





(2.43) 

Notice that all these results give us the eigenvalue +1 of the operator D (see ()2.25|) ). Therefore, in 
one run, if we succeed to perform the above rotations of the atomic states according to (|2.29j) . each 
one followed by the detection of the respective state | fj) or | gj) (J = 1, 2) and | g 3 ) or | / 3 ), we get 
the result of the experiment in favor of quantum mechanics. 

If we perform the test for the state ()2.18|) . following the above procedure, it is easy to see that 
we have the four possible outcomes, 



{Kl, 


1 <?1»(^2, 


92))(K 3 , 


1 /3» - 






(tfl, 


1 9i))(K 2 ,\ 


h))(K 3 , 


1 9s)) «- 






(tfl, 


1 fi))(K 2 ,\ 


f2))(K 3 , 


1 /3» - 




-M-M,-), 


(tfl, 


1 fi))(K 2 ,\ 


92))(K 3 , 


1 93}) <~ 




-M,+M,+). 



(2.44) 

Notice that all these results give us the eigenvalue —1 of the operator D (see (|2.25j) ). as it should be. 
See Fig. 2 where we present a scheme of the GHZ experiment. 



3 ATOMS IN A LAMBDA CONFIGURATION 

Consider a three-level lambda atom (see Fig. 3) interacting with the electromagnetic field inside a 
cavity C. The states of the atom, \a), \b) and |c) are so that the \a) # \c) and \a) \b) transitions are 
in the far off resonance interaction limit. The time evolution operator for the atom- field interaction 
U(t) is given by |2U 

U(r) = -e^ ata |a)(a| + ^(e^ ata + l)|6)(6| + ^(e^ ata -l)|6)(c| +^(e^ ata -l)|c)(6| + ^(e^ ata + l)|c)(c|, 

(3.45) 

where a (a') is the annihilation (creation) operator for the field in cavity C, if = 2g 2 r/ A, g is the 
coupling constant, A = u a — uj b — oj = oj a — uj c — to is the detuning where u a , u b and uj c are the 
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frequency of the upper level and of the two degenerate lower levels respectively and lu is the cavity 
field frequency and r is the atom-field interaction time. For ip = tt, we get 

U{t) = -exp (iirata) \a)(a\ +n+|6)(6| + n_[6><c| +II_|c)(6| + U+\c)(c\, (3.46) 

where 

n + = I( e -^ + l), 

n_ = i( e iWa - 1). (3.47) 

Let us first show how we can get a GHZ state making use of three-level lambda atoms interacting 
with a cavity field prepared in state (|2.1|) . Consider the atom Al in the state \iP)ai = \b\) and a 
cavity C prepared in the state (|2.1jl . We now let atom fly through the cavity C. Taking into 
account (|3.4fij) the state of the system Al — C evolves to 



1 

72 V 



)Ai-c==- f =(\bi)\0)-\c 1 )\l)). (3.48) 



Consider now another three-level lambda atom Al prepared initially in the state I&2), which are going 
to pass through the cavity. After this second atom has passed through the cavity, the system evolves 
to 

\i(>,+)ai-A2-c = -^(|&i)l&2)|0> + |ci)|c 2 )|l)). (3.49) 

which is a a GHZ state. If we had started with the cavity in the state (|0) — \l))/y/2 we would get 
the GHZ state 

I V, ~)ai-A2-c = h) I 6 2 )|0)- I Cl ) I c 2 )|l)]. (3.50) 

In the case of ()3.49|) . following the prescription of the previous section, we can see that in the 
GHZ test we have only four possibilities which are presented schematically below, where on the left, 
we present the possible sequences of atomic state rotations through Kk and detections of | bj) or | Cj) 
(j = 1,2) and | g 3 ) or | / 3 ) and on the right, we present the sequences of the corresponding states 
\ipx-> i) where k = 1, 2 and 3 which corresponds to the measurement of the eigenvalue of the operator 
D given by (j!HT5j) . 



(^1, 


1 ci))(K 2 , 


\C2))(K 3 ,\ 


93)) <- 






(^1, 


1 c 1 )){K 2 , 


\b2))(K 3 ,\ 








(^1, 


1 h))(K 2 ,\ 




93)) *- 






(^1, 


1 6i»(^a,l 


C2))(K 3 ,\ 


/s»<- 




-)l^ 2 ,+)le-)- 



where now, 



(3.51) 

IV* ±> = -^(|6 Jfc >±|c fc » J (3.52) 

for k = 1,2. Notice that all these results give us the eigenvalue +1 of the operator D (see (J2.25|) ). 
Therefore, in one run, if we succeed to perform the above rotations of the atomic states according to 
(|'2.29p . each one followed by the detection of the respective state | bj) or | Cj) (j = 1,2) and | g 3 ) or 
I f 3 ), we get the result of the experiment in favor of quantum mechanics. 
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If we perform the test for the state (|3.5UJ1 . following the above procedure, it is easy to see that 
we have the four possible outcomes, 





Ci))(K 2 , 


c 2 ))(K 3 , 


1 /s» - 






(K h 


ci))(K 2 , 




1 93)) «- 




+>h/4->h/i+>, 




h)){K 2 , 


&2»(#3, 


1 /s» - 




-M,->he->, 


(Kt, 






1 93)) <- 




-M,+M,+>- 



(3.53) 

Notice that all these results give us the eigenvalue —1 of the operator D (see (|2.25j) ). as it should be. 
See Fig. 4 where we present a scheme of the GHZ experiment. 



4 CONCLUSION 



We have presented two schemes of preparation of GHZ states and realization of the GHZ test via 
cavity QED. In these schemes the cavity is prepared in a superposition of zero and one Fock states 
which are not affected dramatically by decoherence and which are relatively easy to handle. In the 
first scheme we make use of atoms in a cascade configuration and in the second scheme we make use of 
atoms in a three-level lambda configuration. The advantage of using a cascade atomic configuration 
is that the detection process is easier, however, we have to perform many atomic state rotations. 
The advantage of using the lambda atomic configuration is that we do not have to use many atomic 
state rotations, however, the detection of the atomic states are more complicated. 



A Time evolution operator 
A.l Two- level atoms 

Let us consider a two- level atom interacting with a cavity field, where |e) and |/) are the upper and 
lower states respectively, with u e and Uf being the two atomic frequencies associated to these two 
states and u the cavity field frequency (see Fig. 1). The Jaynes-Cummings Hamiltonian, under the 
rotating-wave approximation, is given by 

H = ha^a + huj e \e)(e\ + +tioJ f \f}(f\ + hg[a\e)(f\ + a^\f)(e\), (A.54) 

where a) and a are the creation and annihilation operators respectively for the cavity field, g is the 
coupling constant and we write 

H = H + H u (A.55) 

where we have settled 

H = ha) a + hu e \e)(e\ + +Uu>f\f)(f\, 

Hj = hg(a\e}(f\ + (J\f)(e\). (A.56) 
Lets define the interaction picture 

|^) = e^s>, (A.57) 
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l h-\^ s ) = H\^ s ). 



thj t \^ I )=V I \^ I ), 



ge- iAt a^ 



ge iAt a 



where 

Then, we get 
where 

and 
Considering 

\Mt)) = Ui(t)\Mo)) = Ui(t)\Mo)h 

we have to solve the Schrodinger's equation for the time evolution operator 

in—r- = Villi, 



A = {uj e — ujf) — uj. 



dt 



where 



and 



That is, 



Ui(t) 



ltee(t) U ef (t) 
U fe(t) Uff(t) 



Uj(0) 



. d , . 
% 'dt U 

iju ef (t) 

ifu fe (t) 
. d 

l dt Uff{t) 



1 
1 



ge iAt au ef (t), 
ge* At au ff (t), 
ge- lAt a)u ee (t), 
ge- iAt Ju ef (t), 



which can be solved easily using, for instance, Laplace transformation, and we g 
Uj(t) = 



e % 2 '(cos /it — i-^ sin /it) 

inn^ p^ i ~ t -1 (sir 



—ige u 2 "f;(sin /it) a 



m 1 



-iga'e " 2 "-^(sin/xt) 



e 2 



(cos vt + iir- sin vt) 



where we have defined 



/ A 2 

/l = + — , 



v = \ g 2 a^a + 
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In the large detuning limit (A ^> g) we have 



/ , ♦ A A 
y 4 2 



' g 2 a)a + 



A 2 



A 
2~ 



g 2 aa) 



1 \ 
g a'a 



(A.69) 
(A.70) 



and we get easily 

U d (t) = e - ilp{a]a+l) I e)(e I +e ilpata \ f)(f \, 

where (p = g 2 t/A. If we neglect spontaneous emission, the above time evolution operator can also be 
obtained from the effective Hamiltonian 



H d = h 9 -a^a(\e)(e\-\f}(f\). 



(A.71) 



The subindexes in (jA.70j) and (jA.71|) are related to the atom-field interaction described by them, 
that is, a dispersive interaction. 

In the case we have a resonant interaction of an atom with cavity field (A = in (|A.67jl ). if 
the field is a very intense field we can treat it classically. That is, in the time evolution operator 
(jA.67|) we set A = 0, and we substitute the creation and annihilation field operators according to 
a — > rje 10 and a* — > r]e~ id where r\ and 9 are c-numbers. Then, we have a semiclassical approach of the 
atom-field interaction in which the field is treated classically and the atoms according to quantum 
mechanics. In this case (jA.67|) becomes 



cos(grjt) —ie %e sm(gr)t) 
-ie~ %e sin(gi]t) cos(gi]t) 



(A.72) 



Now we take = n/2. If we choose grjt = tt/4 we have the rotation matrix 



1 

Re tl = — 1= 

2 ' 4 y/2 



and for grjt = —ir/4 we have 



For grjt = 71/ 2 we get 



and for grjt = —tt/2 



RlL —2L 

2 ' 4 



RlL ZL 

2 ' 2 



1 

V2 



1 1 
-1 1 



1 -1 
1 1 



Re. 



1 
-1 

-1 

1 



Now if we take 9 = tc and grjt = 7r/4 we have the rotation matrix 



1 

V2 



1 i 
i 1 



(A.73) 



(A.74) 



(A.75) 



(A.76) 



(A.77) 



Choosing the proper values of 9, g, rj and t we can get the rotation matrix we need to perform the 
rotation of the atomic states we desire in a Ramsey cavity. In all the sections of this article when we 
mention a rotation of atomic states in a Ramsey cavity we are taking into account (|A.72|) . 
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Just as a remark, consider (??) and assume we have an intense field, that is, we can use a 
semiclassical approach. In this case we set a — > rje 10 and a) — » 7]e~ ld where 77 and 9 are c-numbers 
and a) a — > rf. Defining (pa) a — > iff] 2 = (3/2, (??) reads 



and for (3 = n we have 



,sc 



IK 



e if3 / 2 



-1 
1 



Any arbitrary 2x2 unitary matrix M may be decomposed as [?] 



M 



-ip/2 Q 







J/3/2 



COS • 



sin 



sin | cos ^ 







J8/2 



(A.78) 



(A.79) 



(A.80) 



where a, [3, 7 and 5 are real parameters. Therefore, we can use (jA.72|) and (jA.78|) to get a rotation 
matrix we need. 



A. 2 Three-level lambda atoms 

We start with the Hamiltonian of a degenerate three-level lambda atom (see Fig. 3) interacting with 
a field cavity mode 



H = huja) a+fuu a \a)(a\ +hu>b\b) (b\+fiu> c \c) {c\ +ha(gi\a)(b\+g2\a)(c\)+ha j '(gl\b)(a\+g2\c){a\), (A. 81) 

where \a), \b) and |c) are the upper and the two degenerated lower atomic levels respectively, a (a^) 
is the annihilation (creation) field operator and gi and g% are the coupling constants corresponding 
to the transitions \a) # |c) and \a) # \b), respectively. In the interaction picture, 



V 

where 



xt gi a \a)(b\ +e- iAlt g^ \b)(a\ +h 



e iA2t g 2 a \a)(c\ +e~ iA2t g& 



c){a 



(A.82) 



The time evolution operator 



whose initial condition is given by 



Ai ~- 


= u a 




— 00 


A 2 = 


= u a 


- 


— to 




U a a 


Uab 


Uac 


U(t) = 


Uba 


Ubb 


Ubc 




Uaa 


U c b 


Ucc 



U(0) 



1 
1 
1 



should satisfy the Schrodinger equation of motion 



ift 



dU 
~dt 



VU 



ClUba + C,2U ca C\Ubb + C,2U c b ClU b c + C,2U C 
Cl u aa Cl u ab Cl u a 

C2 U aa C2 U ab C2 U a 



"etc 



(A.83) 



(A.84) 



(A.85) 



(A.86) 



12 



where 



Ci 

C2 



ke iAlt gi a, 
he iA2t g 2 a. 



(A.87) 



Observe that this equation may be grouped in three sets of couple differential equations. One for 
u aa , u ba and u ca , i.e., 



. fc du aa 

ih — - — 

dt 

. fe du ba 

ifl — : — 

dt 

du ca 

%h — : — 

dt 



another involving only u ab , u bb and u cb , 



ih 



du, 



ah 



ih 



ih 



dt 

du bb 

dt 
du cb 

dt 



and, finally, one involving u ac , u bc and u cc , 

. fc du ac 
ih — - — 
dt 

., du bc 

ih — : — 

dt 

du cc 

ih — : — 

dt 

Let us take Ai = A 2 = A = — irj 



= (lU ba + (2Uca, 

= Cl u aaj 

= C2 U aai 

= ClUbb + (2U c b, 

= (Uab, 

= (Uab, 

= (lU bc + (2Ucc, 

— (Uac, 

= (Uac- 



( x = he iAlt gi a = ha ie v \ 
( 2 = he iA2t g 2 a = ha 2 e vt . 

Notice that all the 3 system of differential equations above are of the form 



.dx 
~dt 



i— = aie^y + a 2 e r > t z, 



le 



'it 



le 



r,t 



dy 

dt 

dz 

dt 



= a\x 



1 x, 



a 2 x. 



(A.88) 



(A.89) 



(A.90) 



(A.91) 



(A.92) 



We can solve the systems of differential equations using, for instance, Laplace transformation and 
we have for the degenerate case 



U aa (t) = 



e 2 



yZ/z cos a//^ — i — sin y/jit 
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where 



Uab(t) 
U ac (t) 
U b a(t) 



.e 2 

-i — — sin y/fitai, 



.e 2 . ^_ 
-i — — sm J/j,ta 2 , 



-%a\ — — sm J [it, 



u bb (t) = l + a{a v 



a\ai + a 2 a 2 

t 1 
Ubc{t) = a[a 2 — f 1 — 

a[ai + a 2 a 2 



e i 2' < (i^ sin J^t + J^ cos J^1 j - y ;' 

e~^' U— ^= sin v/z/t + cos Jlst] — 1 
V 2 V" / 



«ca(*) = 



. e 2 



-za 2 sm Jut, 

VI 1 



1 



u c b(t) = atlati— f r 

t 1 1 



e 2 U 



sin + cos — 1 



: i— + f 

ttjai + a 2 a 2 



e 1 2 * ( sin + y^cos \fut \ — J~v 



(A.93) 



i± t t 

/i = — + a±a[ + a 2 a 2 , 

A 2 t t 
i/ = — ^- + ajai + a 2 a 2 . 

It is easy to show that for the non-degenerate case, i.e., 



(A.94) 



we obtain 

U aa {t) 

Uab{t) 

U ac (t) 

U ba {t) 
U b b(t) 



= giai, 

a 2 = g 2 a 2 , 



e 2 
y^ 



Jjicos JJbt — i— sin Jjlt 



e 2 

-i — — sm J\itai\, 



e 2 

-i — — sm Jjita 2 , 



e 2 



— ia[ sm Jut, 

t 1 1 
1 + a[a\ 



v\ a\ai + a 2 a 2 



(A.95) 



e l 2 [i— sinJUit + jT\cosJT\tj - JU[ 



u bc (t) = a\a 2 - f ^ 



a\a[ + a 2 a 2 



e 2 \i- 



2JV{ 



sin y^t + cos J~u[t — 1 



««»(*) = 



• t e 2 • /— 
-ia 2 sm Jut, 

JV> 
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Ucb(t) 
U cc (t) 



A \ 
i — — sin ^JT>\t + cos ^Jv[t — 1 

2 V^i * " / 



a\ai + a2«2 

t 1 1 r .a, / A __ \ 

l + a ] 2 a 2 —= 1 1 — e 2 U— sin 0^ + V^cos ^v 2 t\ 



where 



A 2 t t 

/i = — + aiai + a 2 a 2 , 



4 

A 2 



= — + a\ai + 0:2 c4> 
A 2 

^2 = + + « 2 a;2. 



Returning to the degenerate case, in the large detuning limit, we have 



and we get 



If 



A 2 a x a\ + a 2 a\ A #1 + 

A /u « = — 

2 A 2 A 

A 2 + A I 0! | 2 + I g 2 | 2 , 
- T + ^ = y + £ a f a, 



«ao(*) 
Waft(*) 

U b b{t) 

U bc {t) 
Uca(t) 
U cb (t) 



exp ^ 1 ■ — ^ ^ 2 ■ taa^j , 

0, 
0, 

0, 



1 + 





9i 


2 


1 9i 


2 + 


92 


2 



exp i 



9i I 2 + I £2 | 2 



ta^a — 1 



= 



9i I 2 + I #2 I 2 

#2*#1 



«cc(0 = 



I 9i 

1 + 



#2 





92 


2 


1 9i 


2 + 


92 


2 



■ f. lgi I 2 + \92 \\ t N 
exp « ra'a 



exp ( % ^ 1 — ^ 2 ■ ta^a ) — 1 



(. \_9i ? + \92 \\ t ^ 1 
exp 2 ra'a — 1 



, (A.96) 



(A.97) 



(A.98) 



(A.99) 



we finally have 



91 = ge*», 

92 = ge^, 



Uaa(t) 



exp % 



(2gH f 
exp 1 ■ 



(A.100) 



15 



Uab(t) 
U ac {t) 
U b a{t) 

Ubb(t) 

Ubc(t) 
U C a{t) 
Ucbif) 

Ucc{t) 



0. 

0. 
0. 
1 

2 

e i{<Pl-V2) 



( ig 2 t + \ ' 

exp \i — a'al + 1 



.2^ 



t, 



exp i z-^-a'a I — 1 



-i{ipi—ip2) 



exp I? — a'a I — 1 



exp I z— ^-a'a 1+1 



(A.101) 



which agrees with the result obtained in |24j . 

Now, let us considering the classical limit in the three level lambda atom interacting with two 
dephased modes with the same frequency (see (|A.96|1 ). We see that the matrix elements of the 
evolution operator read as: 



Ubb 



U c b 



Ubc 



where we have used 



1 + e* x e x — — [e- iAf / 2 (cos p \ e | 2 +A 2 /4t 
A 

+* i = 
2^2 | e | 2 +A 2 /4 

1 + e ;e 2 -^[e- lA */ 2 (cos ^2 | e | 2 +A 2 /4t 



siW2 | eh | 2 +A 2 /4t) - 1], 



+i- 



A 



: sin 



2 | e | 2 +A 2 /4t) - 1], 



2^2 | e | 2 +A 2 /4 
e 1 e 2Tr ^[e-^ /2 (cos ^2 | e | 2 +A 2 /4t 



A 



= sin x /2 lei 2 +A 2 /4t) - 1]. 
2^/2 | e | 2 +A 2 /4 V 

ei^^-^[e- iAi/2 (cos ^2 | 6 | 2 +A 2 /4t 



A 



2 J2 | e | 2 +A 2 /4 



sin a/2 | e | 2 +A 2 /4t) - 1], 



giai 

92^2 



e 1 = ee* e \ 



e 2 = ee 

In the high detuning limit we obtain the expression 

1 



102 



Ubb 



(A.102) 



(A. 103) 
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u cb = V^(e^-l), 

u bc = I e *( e *-1), (A.104) 
where = 9i — 9 2 and p> = 2 | e | 2 t/A . Choosing p> = rr/2 and = 7r/2 we get 

(A.105) 




Choosing ip — ir and = n/2 we get 




Choosing ip — n and = — 7r/2 we get 



Choosing ip — tt and = 7r we get 



Figure Captions 



(A.106) 



U = i(°_ 1 I). (A.107) 



(A.108) 




Fig. 1- Energy states scheme of a three- level atom where |e) is the upper state with atomic 
frequency cu e , \f) is the intermediate state with atomic frequency Uf, \g) is the lower state with 
atomic frequency u g and uj is the cavity field frequency and A = (cu e — Uf) — uj is the detuning. The 
transition | /) ^| e) is far enough of resonance with the cavity central frequency such that only 
virtual transitions occur between these levels (only these states interact with field in cavity C). In 
addition we assume that the transition | e) #| g) is highly detuned from the cavity frequency so that 
there will be no coupling with the cavity field in C. 



Fig. 2- Set-up for the GHZ experiment. Atom AO is sent to the cavity RO, where it is prepared 
in a coherent superposition, and to cavity C to prepare it in a state which is a superposition of 
zero and one Fock states. Atom Al passes through the Ramsey cavity Rl where it is prepared in 
a coherent superposition, cavity C and through the Ramsey cavity R2. Atom A2 passes through 
the Ramsey cavity R3 where it is prepared in a coherent superposition, cavity C and through the 
Ramsey cavity R4. Atom A3 is prepared in state ^3), and sent to the cavity C. Once the GHZ state 
has been obtained, the GHZ test is performed making use of the Ramsey cavities Kl, K2 and K3 
and detectors Dl, D2 and D3 as described in the text. 

Fig. 3- Energy level scheme of the three-level lambda atom where \a) is the upper state with 
atomic frequency u a , \b) and |c) are the lower states with atomic frequency u b and u c , uj is the cavity 
field frequency and A = u a — u b — u = uj a — uj c — u is the detuning. 
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Fig. 4- Set-up for the GHZ experiment. Atom AO is sent to the cavity RO, where it is prepared 
in a coherent superposition, and to cavity C to prepare it in a state which is a superposition of zero 
and one Fock states. Atom Al is prepared in state |6i) and passes through cavity C. Atom A2 is 
prepared in state \b 2 ) and passes through cavity C. Atom A3 is prepared in state \g 3 ), and sent to 
the cavity C. Once the GHZ state has been obtained, the GHZ test is performed making use of the 
Ramsey cavities Kl, K2 and K3 and detectors Dl, D2 and D3 as described in the text. 
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